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We study the Glauber dynamics of simple spin systems to identify dynamical scenarios which may be of 
relevance for the deconfining phase transition in heavy ion collisions. 



1. Introduction 

Lattice gauge theory investigations of the de- 
confining phase transitions have mainly been lim- 
ited to equilibrium studies (an exception is the 
work by Miller and Ogilvie [1]). The equilibrium 
transition with two massless quarks is likely sec- 
ond order and it becomes a crossover with two 
light quarks (and the heavier strange quark) [2]. 
In nature the finite temperature phase transition 
is governed by temperature driven dynamics. 

Early universe: We have a slow cooling process 
(10-5-10-6 >> 10-23 Mogt likely, the effects 
of the dynamics are negligible and no signals of 
the transition are observable nowadays. 

Heavy ion collisions - Bjorken's [3] standard 
scenario: In the center of mass frame the inci- 
dent nuclei are Lorentz contracted into pancake 
shapes. They pass through each other and leave 
behind a region of hot vacuum. The heating is 
presumably not slow on the relaxation time scale, 
but usually considered as a quench, i.e., an instan- 
taneous process. Subsequently, the cooling is not 
much slower than the scale of 10-^'^ s. 

Quenching is a process in which the tempera- 
ture in the symmetric phase below Tc is raised 
instantaneously to a temperature in the broken 
phase above Tc. Quenching has been much stud- 
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ied in condensed matter physics. One finds that 
the dynamics of long-wave modes groups theories 
into dynamical universality classes described by 
the same equations of motion. 

Correlated domains emerge and grow with time 
in such a way that the correlation function of a 
generic field (j) has the simple scaling form 

g{r,r',t) - (0(r, t)0(r ', t)) (1) 
= f{\r~f'\/L{t)), L(t)^tP 

where the exponent p depends only on the dy- 
namical universality class. Properties of the 
Fourier transforms, called structure functions or 
structure factors, allow to differentiate between 
the scenarios of nucleation and spinodal decom- 
position. For spinodal decomposition clusters 
grow on every length scale and one finds pro- 
nounced maxima in the low-momentum structure 
functions. Nucleation is dominated by the growth 
of the largest clusters and the maxima of the low- 
momentum structure functions are insignificant. 

In real heavy ion collisions, the phase transition 
will probably be in between a slow equilibrium 
process and a quench. Therefore, we focus on the 
dynamics of hysteresis loops. We supplement the 
results using equilibrium and quenching data. 

We are interested in the response of an ensem- 
ble of Polyakov loops to a change of temperature 
from the symmetric to the broken phase. Glauber 
dynamics is our best shot which allows us to study 
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the influence of the speed of the heating or cool- 
ing process. It includes canonical Metropolis and 
heat bath updating procedures, which imitate lo- 
cal fluctuations of nature. The aim of such a 
study is not to make precise quantitative predic- 
tions, but to identify possible dynamical scenar- 
ios. 

A direct lattice QCD simulation is very CPU 
time intensive. Therefore, the 3d 3-state Potts 
model serves in our simulations as an effective 
spin model, which is supposed to be in the same 
universahty class [4]. A magnetic field allows to 
imitate the effect of finite quark masses [5] . 

In our first step, reported here, the 2d q-state 
Potts model serves as a test laboratory. The 
advantage is that some results are analytically 
known. We simulate the 2d 2-, 4-, 5- and 10- 
state Potts models, corresponding to a weak 2nd 
order (a = 0), a "strong" 2nd order (a = 2/3), a 
weak first order (a = 1, Ae; = 0.03) and a strong 
first order (a = l,Ae; = 0.35) phase transition, 
where a is the critical exponent of the specific 
heat and Aei the latent heat per link. We study 
the hysteresis with change of coupling 



A/3 = 



2(/?inax /^min) 

n/3 



(2) 



This is a dynamics which slows down with in- 
creasing lattice size. At least 640 cycles are per- 
formed per lattice size. This gives an ensemble of 
non-equilibrium configurations. The equilibrium 
is recovered for n/3 — > oo. 

In the following results are reported for the in- 
ternal energy E as an example of a bulk quantity, 
stochastic clusters and structure functions. 

2. Results 

For a properly rescaled internal energy hystere- 
sis cycles for our ri/s = I dynamics are shown in 
Fig.l. The openings are finite volume estimators 
of the latent heat. Performing infinite volume ex- 
trapolations to the latent heat, the resulting fits 
are only accurate for the q = 10 strong first order 
phase transitions. Otherwise, the results are too 
high when compared with equilibrium values, re- 
sulting in a 'dynamical' latent heat for the second 
order transition. 
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Figure 1. Energy hysteresis loops. Prom left to 
right for g' = 2, 4, 5 and 10. 
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Figure 2. Equilibrium results for the largest clus- 
ter surface. 

Fortuin and Kasteleyn introduced stochastic 
clusters, which allow to formulate the Potts mod- 
els as cluster systems. For a fixed configuration 
the clusters are build as in the Swendsen-Wang 
algorithm. We have studied the Glauber dynam- 
ics of many cluster properties and report here the 
behavior of the largest cluster surface, which has 
pronounced peaks close to the phase transition. 
Fig. 2 shows decreasing peak heights for increasing 
q. This is expected for a scenario, which moves 
from spinodal decomposition to nucleation. Fig. 3 
suggests than that for our ri/j = 1 dynamics the 
transition becomes always spinodal. 



With 



'go 



are defined by 



(r,t), go I 



the structure functions 



Sik,t) 



1 



90=0 



3 




Figure 3. Dynamical n/3 
cluster surface. 
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We recorded the low-lying momenta ki = 
2'jrL-'^ni with rii = (1,0) and (0, 1), na = (1, 1), 
ng = (2,0) and (0,2), m = (2,1) and (1,2), 
n5 = (2,2). 

Spinodal decomposition is characterized by an 
explosive growth in the low momentum modes, 
while the high momentum modes relax to their 
equilibrium values. Miller and Ogilvie [1] inves- 
tigated structure functions under quenching for 
SU(2) and SU(3) lattice gauge theory. Here, we 
study the structure functions during hysteresis 
loops and under quenching. 

Fig. 4 shows Ski{(3) hysteresis for our n^j = 1 
dynamics and confirms that this dynamics leads 
for all q values to spinodal decomposition. In 
Fig. 5 the time evolution of several structure func- 
tions after quenching is depicted for the q = 2 
model with an external magnetic field, so that 
the transition is a crossover. 

3. Summary and Conclusions 

Energy hysteresis cycles allow to locate the 
temperature of the equilibrium transition. For 
weak first order as well as second order transitions 
the dynamics generates a rather large 'dynamical' 
latent heat. 

From our analysis of the Fortuin-Kasteleyn 
clusters and the structure functions we conclude: 
The dynamical transition shows always spinodal 
decomposition and the difference between 1st and 
2nd order is washed out. 
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Figure 4. Equilibrium (small peaks) and rijs = 1 
dynamics behavior of the Ski iP) structure func- 
tion. 
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Figure 5. Time evolution for g = 2 and h = 0.01 
after quenching (3 from 0.2 to 0.6. 

Signals of the spinodal decomposition may still 
survive when there is no proper phase transition, 
but a fast crossover. If a dynamics generates re- 
gions of misaligned Polyakov loops, one expects 
an enhanced production of low-energy gluons. 
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